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Exploring the accelerated expansion of the universe, we investigate the general- 
ized ghost dark energy (GGDE) model from the statefinder diagnosis analysis in a 
flat FRW universe. First we calculate the cosmological evolution and statefinder 
trajectories for non- interacting case and then extend this work by considering the 
interaction between dark matter and dark energy components. We show that in the 
non-interacting case the phantom line can not be crossed and also he evolutionary 
trajectories of model in s — r plane can not be discriminated. It has been shown that 
the present location of model in s — r plane would be close to observational value 
for negative values of model parameter. In the presence of interaction between dark 
matter and dark energy, the phantom regime is achieved, the accelerated phase of 
expansion occurs sooner compare with non-interacting case. Finally, unlike the non- 
interacting case, the evolutionary trajectories in s — r plane can be discriminated 
in the interacting model. Like non-interacting model, in the interacting case the 
present location of GGDE model is closer to observational value for negative values 
of model parameter. 
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I. INTRODUCTION 

The recent astronomical data from SNe la , WMAP 13, SDSS 3 and X-ray 
experiments show that our universe experiences an accelerated expansion. The above 
observational data strongly suggest that the universe is spatially flat and dominated 
by an exotic component with negative pressure, the so-called dark energy [5]. Dark 
energy scenario has got a lot of attention in modern cosmology both from theoretical and 
observational point of view. Observationally, The result of WMAP experiment shows that 
dark energy occupies about 73% of the energy of our universe, dark matter about 23% and 
usual baryons occupies only about 4% of the total energy of the universe c2. Although the 
nature of dark energy is still un-known, but the ultimate fate of the current universe is 
determined by this mysterious component. Theoretically, the first and simplest model for 
dark energy is Einstein's cosmological constant with constant EoS parameter w\ = — 1 . 
Cosmological constant faces with the fine-tuning and cosmic coincidence problems 
recent years a plenty theoretical models have been proposed to interpret the properties of 
dark energy jl-ol. 

Almost all of theoretical dark energy models need to introduce new degree (s) of freedom or 
modifying general relativity. However it would be better to consider a model of dark energy 
without a need of new degree(s) or new parameter in its theory. Recently the so-called 
QCD ghost dark energy has been proposed to interpret the dark energy without any new 
parameter or new degree of freedom (loj. The Veneziano ghost has been suggested to 
solve the U(l) problem in low energy effective theory of QCD Uj. The ghost field has no 
contribution to the vacuum energy density in the flat Minkowsi spacetime. However, in the 
case of curved spacetime it has a small energy density proportional to Ag^^^if, where Aqcd 
is QCD mass scale and H is the Hubble parameter [12]- This model does not encounter 
with some unwanted problems such as the violation of gauge invariance, unitarity and 
causality 12|. Since the Veneziano ghost field is totally embedded in standard model and 
general relativity, one needs not to introduce any new degree(s) of freedom or to modify 
the Einstein's general relativity. The present value of energy density of dark energy in 
this model is roughly of order Aq(j^Hq, with Aqcd ~ lOOMeV and Hq ~ lO^^^eV which 
is in agreement with observed value (3 x lO^'^eV)^ for energy density of dark energy 3]. 
This numerical coincidence is adsorbent and gets rid the model from fine tuning problem 
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10| . Observationally, the ghost dark energy model has been fitted by astronomical data 



including Snia, BAO, CMB, BEN and Hubble parameter data [IJ]. The cosmological 
evolution of dark energy in the QCD ghost model has been calculated in 1^ and has been 
resulted that the universe begins to accelerate at redshift around z ~ 0.6. Also the squared 
sound speed of the dark energy for this model is negative, indicating an instability of the 
model against perturbational theory [l^. In 15|, the ghost dark energy (GDE,hereafter) 
has been extended in the presence of interaction between dark matter and dark energy in 
non-flat universe. The reconstructed potential and the dynamics of scalar fields according 



16|. The reconstructed modified gravity 



the exultation of GDE have been investigated in 
for GDE which describes the late time accelerated expansion has been studied in 17|]. The 
statefinder diagnostic of GDE has been presented in [l8]. Form the statefinder viewpoint, 
the evolution of GDE model is similar to holographic dark energy model and present value 
of statefinder parameters in this model is in good agreement with observation 18 1. 
In all above studies (i.e. Refs. [10-18]), the energy density of GDE is considered proportional 
to Hubble parameter as pd = aH . However, the^energy density of Veneziano ghost field in 
QCD is generally in the form oi H + 0{H'^) [19]. Although, up to now, only the leading 
term H has been assumed for energy density of GDE, but the sub-leading term can also 
be important in the early evolution of the universe 20|. Including the second term in the 
energy density of GDE results better agreement with observation in comparison with usual 
GDE model 21 1. Like 22!], we call this model as generalized ghost dark energy (GGDE). 
The energy density of GGDE model is written as pd = otH + where a and /3 are the 
constants of the model. It has been shown that the GGDE model can result a de-Sitter 
phase of expansion and also in the presence of interaction between dark matter and dark 
energy this model results the phantom regime of expansion [wd < —1) j22|. The other 
features of GGDE model have been presented in [23]. 

It is well known that in addition of dark energy component which describe the accelerated 
expansion of the universe, there exist another mysterious component in the universe 
so-called dark matter. The dark matter component can interpret the fiat rotation curve 



of spiral galaxies and also the scenario of structure formation of universe [2^. Since 
the nature of these component are un known and they they have different gravitational 
treatment, therefore their evolution usually considered independent of each other. However 
recent observation from galaxy cluster Abell A596 indicates the interaction between 
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these components 25|]. Also the observational data from SNIa and CMB experiments is 



compatible with interacting forms of dark energy models [26|. However the strength of this 
interaction is not clearly identified j27|. From theoretical viewpoint it is also acceptable to 
consider the interaction between dark matter and dark energy. In the unified models of 
field theory dark matter and dark energy can be interpreted by a single scalar field in a 
minimally interaction. Also considering interaction between dark matter and dark energy 
can solve the coincidence problem [s^ 



In this work our main task is to investigate the interacting GGDE model in statefinder 
diagnostic analysis. Different dynamical dark energy models obtain accelerated expansion 
at the present time (g < 0), where q is deceleration parameter. Hence we need a diagnostic 



tool for discriminating these dark energy models. For this aim, Sahni et al. [28[ and Alam 
et al. {29!, by using the third time derivative of scale factor, introduced the statefinder pair 
{s,r}. These parameters in flat universe are given by 



'd r — 1 

3(g-l/2) 

The parameters s and r are geometrical, because they only depend on the scale fac- 
tor. In statefinder analysis we plot the evolutionary trajectories of dark energy model is 
s — r plane. In recent years the various dark energy models such as quintessence, holo- 
graphic,new holographic, phantom, tachyon, chaplygin gas, agegraphic, new agegraphic, , 
polytropic gas and ghost dark energy models have been studied in the statefinder analy- 



sis [28, 



29 



31 



32| . These models have different evolutionary trajectories in {s, r} plane. 



therefore the statefinder tool can discriminate these models. The standard ACDM has no 



evolution in this plane and corresponds to the fixed point {s=0,r=l} 28|. The present ob- 
servational value for statefinder parameters are {sq = — 0.006, r = 1.02} j3o|. The distance 
of the current value of statefinder pair {sq^tq] of a given dark energy model from the ob- 
servational value {sq = — 0.006, r = 1.02} is a valuable criterion to examine of model. Here 
we see that the location of standard ACDM model in s — r plane is near to observational 



value. In 



321], the evolution of original GDE has been calculated by statefinder diagnostic 



in s — r plane and shown that the GDE model mimics the ACDM at the late time. Also the 
behavior of GDE is similar to holographic dark energy in this plane 32|. 



In this work we first calculate the cosmological evolution of GGDE model and then investi- 
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gate this model from statefinder diagnostic analysis. The paper is organized as follows: In 
sect. 2, The GGDE model is presented in non-interacting universe. The interacting case of 
GGDE model is given in sect. 3. We calculate the numerical results in sect. 4 and conclude 
in sect. 5. 



II. NON-INTERACTING GGDE MODEL 

A flat Friedmann-Robertson- Walker (FRW) universe dominated by dark matter and dark 
energy is given by 

= + (2) 

where pm and pd are, respectively, the energy density of pressureless dark matter and dark 



energy and nip is the reduced planck mass. The energy density of GGDE is given by 22| 



pj^ = aH + I3H^ (3) 

where a and /3 are constants of model. The Friedmann equation ([2]) in terms of dimensionless 
parameters is written as 

^m + ^A = l- (4) 

where 

Q _ Pm _ Pm Q _ P£ _ Prf /_\ 

The conservation equations for pressureless dark matter and dark energy without interaction 
read the following equations 

Pm + SHpm = 0, (6) 

Pd + 3/7(1 + Wd)pd = 0. (7) 
Taking the time derivative of Friedmann equation ([2D and using (jH [6], [7]) obtains 

§-, = -l[^ + y^A^d] (8) 

Differentiating Eq.([n]) with respect to time yields 



p = H{a + 2f3H) 



(9) 
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Inserting (|9]) and ([3]) in conservation equation for dark energy ([7]) and using ([8]), the EoS 
parameter of GGDE model can be obtained as 

where ^ = SirG^/S. In the limiting case ^ = 0, this relation reduces to its original form in 



15|. The deceleration parameter q by using and (ITUI) . in GGDE universe is obtained as 



_ 1 3 ^-Qd . . 

^ 2 2{^ + Q,-2) ^ ' 

The decelerated phase of expansion at the early time is indicated by g < and accelerated 
phase is related to g > 0. Taking the time derivative of dimensionless dark energy density in 
(|5]) and using ([3]), iQ, we obtain the equation of motion for the evolution of energy density 
of GGDE model as 

where prime is derivative with respect to In a. Taking a derivative of ( ITOj) with respect to 
In a, the equation of motion for EoS parameter can be calculated as 

, 3(l-f],)(e-fi,)r , (^-^]^)(2-2^],-0- 
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(13) 



Using the above relation, in this stage, we calculate the statefinder parameters s and r for 
GGDE model in non-interacting universe. In general form, relation ([T]) for a given dark 
energy model in flat universe can be written as 

9 3 
r = 1 + -Wd^ail + Wd^d) - - (w'^^d + Wd^'a) (14) 



and 



Inserting relations (1121) and ( 1131) in equations ( IHI) and ( jlSl) . we obtain the statefinder pa- 
rameters for GGDE model in spatially fiat universe 



s = i + Wdnd-]-i^ + ^) (15) 

3 Wd iid 



(16) 
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In the limiting case of dark energy dominated universe {Qd ^ 0) the parameters {s,r} 
tends to {0, 1}, respectively. Hence the GGDE model mimics the ACDM model at the late 
time when Qd 0. In sect. IV, we calculate numerically the evolution of GGDE model in 
non- interacting universe from the statefinder viewpoint. 



III. INTERACTING GGDE MODEL 



In this section we consider the interaction between dark matter and dark energy compo- 
nents. In this case the conservation equations for these components are: 

p'm + SHpm = Q, (18) 
Pd + 3H{1 + Wd)pd = -Q. (19) 

where Q in right hand side indicate the interaction term. The positive value of Q means the 
transition of energy from dark energy to dark matter component. It should be noted that 
the left side of ([6]) and (JTj) are inversely proportional to time. Therefore the parameter Q 
can be considered as a function of Hubble parameter H such as following forms: 

(i) Q oc Hpd 

(ii) Q oc Hpm 

(iii) Q oc H{pm + pd)- 

One can assume the above three forms as Q = Fp^, where for case (i) F = 3b^H, for case 
(ii) F = Sb'^Hjf^ and for case (iii) F = Sb'^H-^. The parameter 6 is a coupling^ constant 



indicating the strength of interaction between dark matter and dark energy [33[. In this 
work we assume the third form of interaction (i.e., Q = SHb"^-^). 

Substituting Q in f|T9|) and using ([3]), ([8]) and Qj, the EoS parameter of interacting GGDE 
model is obtained as 

e - fid - 262 

Inserting 6 = recovers the EoS parameter of non-interacting case in previous section. 
Substituting (l20l) in ([8]) results the deceleration parameter q for interacting follows 

3 {^~Qd-2b') 
^ 2 2 (e + fid-2) ^ ^ 

It has been shown that for selected parameters (^ = 0.03, b = 0.15, Qdo = 0.72) the deceler- 



ation parameter at the present time is go = —0.38 which is consistent with observation 



22|. 
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Taking the time derivative of dimensionless dark energy density in ([5]) and using ([3]), ([9]) 
and fl20|) , the equation of motion for the evolution of energy density of interacting GGDE 
model can be obtained as 



n'=-3 



(22) 



In the limiting non-interacting case (6 = 0) the respective relation in previous section is 
retained. Derivative of ( 120|) with respect to In a results 



3(1 fid -262) + 62(2 _f]^_^)r (^_^^_262)(2-2fi,-0 



1 + 



fi,(2 - firf - 



(23) 



fi,(2 -ild- 0' 

Once again, the respective relation in previous section can be obtained by setting 6 = 0. 
Finally by substituting relations fl22|) and fl23|) in general relations f|T^ and f|T5|) . we obtain 
the statefinder parameters for interacting GGDE model as follows 

fe-fid-262)(l-fi,-62) 



1 + 9- 



(2 - fid - 0' 



:i - fin 



- fid)(e - fid -26^) + 6^(2- fid -0 ^^ , ^ 



(2-fid-O' 



:i - fid) + 6' 



,2 2-f2d-g 



e + 6^ 



(24) 



(25) 



(2-fi,-0 (2-fi,-02 ' 

Setting 6 = retains the relations for s and r in previous section. In next section we 

investigate the evolution of interacting GGDE model in s — r plane can calculate the effect 

of interaction parameter 6^ on the evolution of the model. 



IV. NUMERICAL RESULTS 



Here we present numerical description for cosmological evolution and statefinder analysis 
of GGDE model in the fiat FRW cosmology. In numerical procedure we fix the cosmological 
parameters at the present time as fi^ = 0.3 and = 0.7. We first consider non-interacting 
case and then interacting case of GGDE model. 



A. non-interacting case 



The EoS parameter of non-interacting GGDE model as a function of density parameter 
fid is given by f lTO|) . By solving coupled equations f|T2|) and f lTO|) . the evolution of EoS 
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parameter in terms of cosmic redshift z = 1/a — 1 and for different illustrative values 
of ^ is shown in Fig.(l). The cosmic redshift z = represents the present time, z > 
indicates the past times and z < expresses the future. We see that for any value of ^ the 
non-interacting GGDE model can not enter the phantom regime {wd < —1) at all. We also 
see that the EoS of GGDE model with ^ > is larger than EoS of standard GDE model 
= 0.0). At the future epoch, the EoS tends to —1 which implies that the GGDE model 
mimics the cosmological constant at that time. 

The deceleration parameter q which indicates the decelerated or accelerated phase of 
expansion for non-interacting case is given by ( fTTl) . Solving coupled equations ( JTTl) and 
( [T2|l . the evolution of parameter g as a function of redshift parameter z has been shown in 
Fig. (2) for different values of model parameter ^. The standard GDE model is indicated by 
solid line, one can see that for ^ < the GGDE model enters the accelerated phase sooner 
and for ^ > later compare with standard GDE model. 

The statefinder pair {s,r} for non-interacting GGDE model is given by relations ( |T6ll 
and f|T7|) . Solving these coupled equations together with f|T2|) . we obtain the evolution of 
parameters s and r in terms of redshift z. In Fig. (3), we plot the evolutionary trajectories 
of non-interacting GGDE model for different values of model parameter ^ in s — r plane. 
We see that the parameter r first decreases and then increases and also the parameter s 
decreases during the history of the universe from past to future. The important note is that 
in non-interacting case the GGDE model has been shown by single evolutionary trajectory 
for any value of ^. Hence the evolutionary trajectories can not discriminated by model 
parameter ^. The present value of statefinder pair {sq, vq} is indicated by colored circle on 
the figure. Also the location of ACDM model in s — r plane ,i.e., (s = 0,r = 1), has been 
shown by star symbol. The other feature is that the present value {scTq} is discriminated 
by model parameter ^. In the case of < 0, the distance of {sq, tq} from observational point 
{so = — 0.006, r = 1.02} (red star point on the figure) is shorter compare with standard 
GDE model (i.e., ^ = 0.0). While for ^ > the distance is larger than GDE model. 



10 



B. interacting case 

Here we calculate the numerical description of cosmological evolution and statefinder 
diagnosis for interacting GGDE model in spatially flat universe. First the evolution of EoS 
parameter in terms of cosmic redshift is plotted in Fig. (4). For this aim we solved the 
coupled equations ( 120|) and (1221) . In left panel, by fixing interaction parameter as 6 = 0.2, 
the EoS parameter Wd is plotted for different illustrative values of model parameter as 
described in legend. The solid line represents the original GHDE model. For all cases of 
^, the interacting case of GGDE model can cross the phantom line {wd = —1) form upper 
limit(wd < —1) to lower limit {wd < —1). This behavior of interacting GGDE model in 
which the phantom line is crossed from up to below is in agrement with recent observations 



351]. In right panel the model parameter is fixed as ^ = 0.1 and the interaction parameter b 
is varied. We see that the non-interacting GGDE model (=0.0) cannot enter the phantom 
regime (solid line), while in other cases ^ 0) the phantom regime has been achieved. 



In Fig. (5), by numerical solving of relations ( l22l) and ( 12T|) . the evolution of deceleration 
parameter in terms of redshift has been shown in the context of interacting GGDE 
model. In left panel, by fixing b = 0.2, the parameter ^ is varied as indicated in legend. 
Same as non-interacting case, the parameter q starts from positive value at the earlier 
( representing the decelerated phase at the past time) and ends to negative value later 
(indicating the accelerated phase at the present time). Like non- interacting case, transition 
from decelerated phase to accelerated phase for < takes place earlier and for ^ > 
later, compare with original GDE model (see left panel of Fig. (5)). In right panel, for an 
illustrative value ^ = 0.1, the evolution of q has been shown for different values of interaction 
parameter b. The interaction parameter b can influence on the transition epoch from g > 
to g < 0. We see that in the framework of interacting GGDE model the accelerated phase of 
expansion {q < 0) can be achieved sooner for larger values of b. For all cases, g — > — 1 at the 
late time indicating that the GGDE model mimics the standard A CDM model at that time. 

Finally, by solving relations (125|) and [2^ . the evolution of interacting GGDE model in 
s — r plane is plotted in Fig. (6). In left panel, by fixing ^ = 0.1, the interaction parameter b 
is varied as indicated in legend. The evolutionary trajectories starts from right to left and 



11 



ended at the ACDM fixed point. The distance of present value {sq, tq} form the observational 
point {sq = —0.006, r = 1.02} (red star point on the figure) becomes larger by increasing b. 
In right panel, by fixing b = 0.2, the trajectories have been plotted for different illustrative 
values of ^. The important note is that, contrary with non-interacting case, in the presence 
of interaction between dark matter and dark energy (6 ^ 0) the evolutionary trajectories in 
s — r plane are discriminated by parameter ^. Also the distance of {sq, tq} from observational 
point {sq = —0.006, r = 1.02} is shorter for ,^ < and larger for > compare with original 
GDE model (i.e., ^ = 0.0). 

V. CONCLUSION 

In summary, we considered the generalized version of QCD ghost dark energy (GGDE) 
model in both non-interacting and interacting universe. The cosmological evolution and 
also the statefinder diagnosis of the model have been calculated. We showed that: 

(i) . In the non-interacting GGDE model the phantom regime can not be achieved and the 
EoS parameter reaches to asymptotic value = — 1 at the late time. We also showed that 
for negative values of model parameter ^ the transition from decelerated to accelerated 
phase takes place sooner compare with original ghost dark energy (GDE) model. The 
statefider analysis was also performed for non-interacting GGDE model. In this case, in the 
absence of interaction between dark matter and dark energy, the evolutionary trajectories 
of model in s — r plane can not discriminated. However, the present value of statefinder 
parameter {so,ro} of the model is diagnosed by parameter ^. We concluded that {so,ro} is 
closer to observational value {sq = — 0.006, r = 1.02} for negative values of ^ (see Fig.(3)). 

(ii) . In the presence of interaction between dark matter and dark energy, the GGDE model 
can cross the phantom line from up to down in agreements with observations [35]. Also, 
in the context of interacting GGDE model the esterase to accelerated phase of expansion 
occurs earlier compare with non-interacting case. Finally the statefinder diagnosis analysis 
was also performed for interacting case and we showed that in the presence of interaction the 
evolutionary trajectories of GGDE model s — r plane are diagnosed. Same as non- interacting 
model, the present value {scTq} is closer to observational point for negative values of ^ 
(see Fig. (6)). 
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FIG. 1: The evolution of EoS parameter of non- interacting GGDE model in terms of redshift 
parameter z for different illustrative values of model parameters ^. Here we take = 0.70 and 
Vt, = 0.30. 




FIG. 2: The evolution of deceleration parameter q in the context of non-interacting GGDE model 
as a function of redshift z for different illustrative values of model parameters ^. We toke 0|J = 0.70 
and = 0.30. 
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FIG. 3: The statefinder plot of non- interacting GGDE model in flat FRW universe with = 0.70 
and $7^ = 0.30. The present value {sq, vq} is indicated by colored circle on the curves. The location 
of ACDM model and observational point are indicated by black and red stars, respectively. 
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FIG. 4: The evolution of EoS parameter of interacting GGDE model versus redshift parameter z 
as described in legend. In left panel, the interaction parameter b is fixed and in the right panel the 
parameter ^. 
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FIG. 5: The evolution of parameter q for interacting case of GGDE model as afunction of redshift 
z as indicated in legend. 
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FIG. 6: The statefinder plot of interacting GGDE model. In left panel the parameter ^ is fixed 
and in right panel the parameter b is fixed. The present value {sq, tq} is indicated by colored circle 
on the curves. The location of ACDM model and observational point are indicated by black and 
red stars, respectively. 
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